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In satellite systems in which the relative positions of satellites are al- 
lowed to vara, there will he periods during which no service will be pro- 
vided between given ground stations. Such periods are called "outages," and 
the intervals between successive outages are called "innages." Here the 
outage and innage time distributions are studied with the help of an anal- 
og!/ between a satellite system and a traffic system. The arrival of a cus- 
tomer in the traffic system corresponds to a satellite coming into view, and 
the service time of the customer corresponds to the time the satellite remains 
in view. In particular, the methods of analysis developed for traffic systems 
are applied to determine an approximation for the distribution of innage 
lengths. 

I. IXTHODUCTION 

Several types of communication systems have been proposed which 
would use repeaters orbiting the earth as artificial satellites. The prob- 
lem considered in this paper arises in systems employing a number of 
satellites at altitudes of several thousand miles. Typically, the orbit 
altitude might be of the order of ";()00 miles with a period of revolution 
of about ") hours. 

The companion paper by Rinehart and Robbins 1 discusses the condi- 
tions under which a particular satellite will be visible to a given pair of 
ground stations. For the orbit altitudes considered here, the satellite 
will be visible intermittently. Conceivably the relative positions of the 
satellites could be maintained so that at least one satellite is mutually 
visible from the two ground stations at all times. However, at least for 
some of the early systems proposed, it is of interest to consider the case 

1(571 
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in which small differences in orbital period cause the relative positions 
of the satellites to vary with time. Attention is therefore directed to the 
statistical characteristics of satellite visibility. 

We are especially interested in those periods during which no satellite 
is available for communication between a given pair of ground stations. 
For convenience these events are called outages although, as pointed out 
by Rinehart and Robbins, these occurrences need not imply an inter- 
ruption of calls in progress. By analogy the intervals between outages 
are called innages. 

For any particular system the most effective method of obtaining 
outage information appears to be that of simulation. The course of the 
system for a year or more is computed, with the help of a high-speed 
digital computer, and the outages and innages recorded. [A method for 
doing this and some sample results are presented in the paper by Rine- 
hart and Robbins]. 

The present paper is concerned with determining the distribution of 
outage and innage lengths. The theory of outage length distribution has 
been studied by a number of people and some of their results, namely 
those which are needed here, are summarized in Section II. The traffic 
model is described in Section III. Sections IV and V contain results 
predicted by the model. In Section VI the predicted results are com- 
pared with those obtained by simulation. The work of Appendix B gives 
the basis for a model which is simpler but less accurate than the one 
described in Section III. The material in both appendices also appears 
to be of interest from the standpoint of traffic theory. 

The general conclusion is that, for the cases examined, both the out- 
age and innage distributions are approximately exponential. Their aver- 
ages are related by the rather simple equation (1). 

II. PRELIMINARY RESULTS 

First note that the average innage length h and the average outage 
length t ("i" for innage and "o" for outage) are connected by the re- 
lation 

h = r^-t (l) 

1 - q 

where q is the quality of service; i.e., the fraction of time transmission 
is possible. This relation is always true and follows almost immediately 
from the definition of q. Typical values are q = 0.99 and t = 0.25 hours; 
and it follows that the corresponding average innage length is U = 
24.75 hours. 
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Next, let k be the number of satellites in the system. Consider a 
particular satellite. Every now and then it will pass over the region of 
mutual visibility. Let 

T = average length of time between its reappearances. For polar 
orbits, T is the orbit time if the region of visibility includes 
one of the poles. 

b~ = average length of time the satellite is visible during one pass. 
The quantity 6 is a rate which occurs frequently in the analysis. 

p = fraction of time the satellite is visible. From these definitions 
it follows that 

6" 1 = pT. (2) 

The satellite systems considered here are restricted to those for which 
T and b~ are almost the same for all k satellites. 
With this notation we have 

g-l-(l-p)*, (3) 

1= (1 - p) T/k, (4) 

where P(t) is the probability that (he length of an outage will exceed t. 
Also the expected number of satellites visible at a given time is kp. 
These formulas are based upon the assumption that the phase angles of 
the satellites are independent random variables and are distributed uni- 
formly over the interval (0,27r). Equation (3) is due to J. R. Pierce and 
R. Kompfner." Equations (4) and (5) have been given by R. E. Mosher 
and R. I. Wilkinson, respectively, in unpublished memoranda. For 
values of / and k of interest, (5) may be approximated by 

P(0 = exp (- t/t ). (6) 

It will be observed that the constants of the satellite system enter the 
right-hand sides of (3) to (5) only through the three parameters k,p,T. 
Thus, as far as these formulas are concerned, the satellite system is 
specified by k,p,T. 

III. TRAFFIC SYSTEM MODEL 

The satellite system will be represented by a traffic system model 
which consists of k independent servers, each having an average service 
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time 6 -1 and a service time distribution function B{t). Customers are 
supplied to the servers by k independent Poisson sources, each producing 
customers at an average rate a. When a source produces a customer, that 
particular source is removed from the system while the customer is 
being served. Thus, when n servers are busy serving n customers, the 
average arrival rate is (k — n) a. This type of input is a special case of a 
more general type (the limited source or "Engset input") which has 
been studied in traffic theory. 

The instant a satellite becomes visible from both the receiver and 
transmitter corresponds to the arrival of a customer. The length of 
time a satellite remains visible corresponds to the time required to serve 
a customer. After a customer has been served he leaves the system. This 
corresponds to the satellite leaving the region of mutual visibility. The 
state in which n satellites are visible simultaneously (state n) corre- 
sponds to the state in which n servers are busy serving n customers. An 
outage corresponds to an idle period (state 0), i.e., a period during which 
all servers are idle. An innage corresponds to a busy period, i.e., to a 
period when one or more servers are busy. 

Note that the constant orbit time of the satellite may introduce a 
regularity in the satellite arrivals. The traffic model has the shortcoming 
that there is no corresponding regularity in the customer arrivals. 

The analogy between the satellite system and the model is established 
by taking k and b~~ l to have the same values in both and setting 

a = i/(i - p )T. (7) 

To justify this choice for a, note that if a particular satellite is not visible 
at a time t selected at random, the chance that it will become visible in 
t, t + dt is dt/(T — pT). Comparison with the corresponding proba- 
bility a dt for the traffic system gives (7). When the three satellite sys- 
tem parameters k,p,T are known, the three model parameters k, b~ , a 
follow at once from (7) and fe _1 = pT. 

IV. OUTAGE AND QUALITY OF SERVICE PREDICTED BY MODEL 

The values of the quality of service q and average outage length t 
predicted by the model agree exactly with (3) and (4) while the pre- 
dicted outage length distribution is the exponential approximation (6) 
to Wilkinson's polynomial expression. A sketch of the proof of these 
statements is given in the following paragraphs. 

First assume the service distribution B(t) to be exponential, i.e., 
equal to 1 — e' bl . Then since the sources are Poisson, the behavior of 
the system is governed by the k + 1 state equations (Ref. 3, p. 30) 
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P ' = - fettPo + bl\ 

1\' = kaP - [(k - l)a + b]Pt + 2bP 2 

P 2 ' = (k - l) a p 1 - [(k - 2) a + 2b]P 2 + 36P 3 
P,' - (A: - 2)aP 2 - [(A: - 3)a + 36]P 8 + 46P 4 



(8) 



P K ' = aPk-x - kbP k 

whore P„ bs P„(/) is the probability the system is in state n at time t 
and primes denote time derivatives. The steady state probability p n 
that, exactly n customers are present at a time picked at random is 

P " = Ci) (a/rl) " 7 '"' Po = (1 + ab" l )- k . (9) 

This follows upon setting the derivatives in (8) to zero, taking P„(t) = 
p H , and solving the equations step by step. The expected number n of 
customers present is kab~ l /(l + a6 _1 ) and the average arrival rate is 

a(k - n) = /••«/(! + ab' 1 ). (10) 

The (juality of service is now 

9=1- Pa (11) 

and expression (3) for q may be obtained by using a = 1/(1 — p)T, 
b~ l = pT. Since, (i) an outage corresponds to state 0, («') state can 
end only through an arrival, and (Hi) the arrivals in state are Poisson 
with rate ka, it follows that exp (— hat) is the probability (predicted by 
the model) that the length of an outage will exceed /. This agrees with 
the exponential approximation ((i), and the average value ~t a = l/(ka) 
agrees with (4). 

It is known that expression (9) for the steady state probability 
p n holds not only for exponential service but also for the general service 
distribution B(t) (Ref. .3., p. 90). Hence the model predicts that expres- 
sions (3), (4) and (6) still hold when the length of time a particular 
satellite stays in view has an arbitrary distribution B(t). 

V. INNAQE DISTRIBUTION PREDICTED BY MODEL 

The average innage length (i.e., average busy period) predicted by 
the model when the service distribution B(t) is arbitrary follows from 
t = l/ka and q = 1 — p : 
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U = X -^^l = [(1 + ab~ x ) k - l]/ka. 

This much is easy. It is much more difficult to get the complete distribu- 
tion, as the following work shows. 

When the service lengths are exponentially distributed, the busy 
period distribution may be obtained by solving a "first passage" prob- 
lem. State is made absorbing and the system is started with an arrival 
at time 0. Thus the system starts in state 1 at time and jumps from 
state to state in accordance with the arrivals and departures of custom- 
ers. The system eventually lands in state and stays there. This cor- 
responds to the end of the busy period or innage. 

When state is made absorbing, the first two of the k + 1 state equa- 
tions (8) are replaced by 

( 12 ) 
TV = _ [(ft - i)« + b]P! + 2bP,. 

The modified equations (8) are to be solved subject to Pi = 1 and 
Po , P2 , ■ ■ • , Pk = at time 0. The probability that the length of an 
innage will exceed the length t is 

(7(0 = 1 -Po. (13) 

Step-by-step computation of the derivatives of Pa at t = from the 
differential equations gives the power series 

\(k - l) 2 a + 4( A: - 1) ab + tf] b? 
3!~ (14) 

+ [(/>• - 1)V+ (k - l)(9/o - 11) ab 

+ 1KA-- l)«6 2 + b 3 ]^+ ••• 

which is useful for small values of t. 

Since Pi is determined by the last k differential equations of the 
modified set, and since the coefficients in these equations are constants, 
we may expect Pi [and hence G(t)] to be expressible as the sum of k 
exponential terms. Indeed, when ^„(s) is used to represent the Laplace 
transform of P„(0» the k differential equations for Pi(0, ■ • ■ , P*(0 go 
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into k linear equations for <pi(s) , ■ ■ ■ , <p k (s). By introducing the genera- 
ting function 

<f>(.r, s) = .(-^(s) + • • • + x k Vk (s) 

in the usual way, the linear equations may be combined into 

Ms) - * = (1 - a;) (6 + ax) 3 *d*' S) - (s + ak - akxMx, s). 

To obtain <pi(s), rewrite this equation as 

[Ms) - -f](l - *r*(6 + ax)—*" 1 

5 , (15) 

= f-(l - *)'<& + ax)-'~ **(*, s) 
ox 

where z = s/(a + b). Assume Re (2) > and integrate (15) from x = 
to x = 1. The right-hand side vanishes because <£(0, s) is zero. Changing 
the variable of integration on the left-hand side from x to y = (1 — x)/ 
(b -+- ax) gives 

Ms) f ' /T'U + ay) h dy = f (1 - s/W'O + otyf^dy. (16) 

Expanding (1 + ay) k and integrating termwise shows that the coefficient 
of 6<£i(s) is b~*Fk(z) where 

F t (z) = ± (';) ( -^T . (17) 



„=o \n/ z + 1 

When the integrand on the right-hand side of (10) is replaced by its 
equivalent 



[ z(a + b)1 --i M . v* a + /> f/ ... . 



(which is suggested by an integration by parts and adding and sub- 
tracting various terms) we obtain 

bp 1 (s)F k (z) = (l + ±) F k (z) - A (1 + oft" 1 )* 1 . 

From P ' = bPi it follows that s<p t) (s) = b<pi(s) and hence 

1,1 6(1 + ah- 1 ) 1 * 1 



<Po(s) = - + r 

S /.-a kasPk(z) 

As s -> oc, *•*(*) tends to (1 + a/;-')*'/* and <p„(s) is 0(l/s). Writing 
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<p (s) as the sum of partial fractions, inverting to obtain P (t), and using 
(13) gives the expression we seek: 

„, , _ « 6(1 + ab- l ) k+l exp[(q + b)z m t\ . . 

W) ~ £„ kaz m F k '(z m ) 

where Z'YO) = dF k (z)/dz and z , z x , • ■ ■ , z k -x are the zeros of F k (z). 
These zeros lie between the poles at 0, -1, -2, • • • , -k. When t is 
large G(t) is given, effectively, by the term corresponding to m = in 
(18). Usually zo is close to 2 = and may be obtained by successive 
approximations from F k (z) = 0, i.e., from 

1 £ /k\ (ab- i r 

The foregoing method is the one originally used to obtain G(t) as the 
sum of exponential terms. Subsequently, a more elegant method of ob- 
taining (18) for (7(/) was suggested by L. Takacs. His method removes 
the restriction that the service time distribution be exponential. 

Takacs' result is the following: Let /3(s), y(s) be the Laplace trans- 
forms of B'(t), —G'(t), the service time and busy period probability 
densities, respectively. The equation to determine y(s), given /3(s), is 



= f e- s '[P att (()} k dt 



19 



s + ka — kay(s) 

where P m (t) has the Laplace transform l/[s + a - a(3(s)]. It will appear 
later that the subscript a refers to the idle state (state 0). The expres- 
sion (18) for G(t) may be obtained from (19) by starting with /3(s) = 
b/(s + b). It turns out that P aa (t) is given by [b + a e ( bu ]/(b + a) 
so that the Laplace transform of [P„„(0f is not difficult to compute. 

One way to establish (19) is to regard the model as composed of k 
independent simple systems, each consisting of a source connected to a 
server. Consider a simple system. The lengths of the idle and busy 
periods have the respective probability densities a e~ a , B'(l) with 
Laplace transforms a/ (a -\- s), (}(&). 

From (26) of Appendix A, the probability P na {t) that an idle period 
will be in progress at time t, given that one is in progress at time 0, may 
be determined by inverting its Laplace transform l/[s + a — a8(s)]. 
The probability that all A- servers are idle at time t, given that they are 
idle at time 0, is [P a „(of- Equation (19) now follows upon using (26) 
again. This time the type (a) intervals correspond to the periods (out- 
ages) during which all k servers are idle. The arrival rate is ak and p a (t) 
is ak e~ akt . The type (b) intervals correspond to the innages with 
probability density -G'(t) and Laplace transform y(s). 
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VI. COMPARISON WITH SIMULATION" 

It is interesting to compare the results predicted by the model with 
those obtained by simulation. As an example, we shall take results ob- 
tained by Rinehart and Robbins for a Maine- Western Europe link. This 
link was assumed to have 18 satellites in random polar orbits at a 
height of 0,000 nautical miles. 

The positions of the orbits, the locations of the receiver and trans- 
mitter, and some computations involving a number of representative 
passages over the region of mutual visibility lead to 

1. the value T = 6.35 hours for the average interval between reap- 
pearances of a particular satellite, 

2. the distribution function B(t) for /, the length of time it is visible, 

3. an average value of / equal to b~ l = 1.46 hours, 

4. the value p = b~ x /T = 0.230 for the fraction of time the satellite 
is visible. 

It turns out that the probability density B'(<) can be approximated by 
the rectangle 

, t < 0.86 hours 

B'(t) = - 1 0.833, 0.86 < / < 2.06 (20) 

, 2.06 < /. 

The system was required to furnish a quality of service close to 0.00. 
This requirement together with the expression q = 1 — (1 — p) k and 
the value p = 0.230 gives A* = 18 and q = 0.09094. The model param- 
eters are therefore taken to be 

fe = 18, /;"' = 1.46 hours, a ' 



(20 

= = 205 

(0.770) (6.35) 

with the understanding that p = 0.230. 

The values of q, t„ ,h obtained by simulation are compared with those 
predicted by the formulas of Section II (and also by the model) in 
Table I. 

The second column gives values obtained by Rinehart and Robbins 
by a simulation which followed the system for 18 months. The third 
column gives values computed from the model parameters (21). It is 
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Table I 



Quantity 



U 

b ' 



Simulation 



0.9918 
0.291 hours 
35.1 hours 



Model, p = 0.230 



0.99094 

l/ka = 0.271 

q(l-q)- r t u = 29.6 

1 .46 hours 



Model, p = 0.234 



0.9918 
0.270 
32.6 
1.48 



seen that the values of the average outage length t agree fairly well. 
The discrepancy in U reflects the shortcomings of the model. The last 
column shows what happens when we hold k and T at their former values 
of 18 and 6.35, and fudge the value of p so as to make q have the simula- 
tion value 0.9918. This changes p from 0.230 to 0.234 and U from 29.6 
to 32.6. It is seen that the value of h is very sensitive to such changes. 
Fig. 1 shows three curves for the innage length distribution. The ordi- 
nate is G(t), the probability that an innage length will exceed t. Curve A 
is the curve predicted by the model assuming exponential service, and 
is obtained by substituting the parameter values (21 ) in expressions ( 14) 
and (18) for G(t). Curve B is the exponential approximation 

0(0 - <r l,li (22) 

with h equal to the simulation value 35.1 hours. Curve C is the result 
obtained by simulation. During the 18 months simulated there were 122 
innages, the longest lying between 250 and 260 hours, the next longest 
between 190 and 200 hours, and so on. The average innage lengths cor- 
responding to curves A and C are h = 29.6 and U = 35.1, respectively, 
in agreement with the table above. 

The curves shown in Fig. 1 agree moderately well. The agreement 
would be improved if curve A could be shifted so as to give an average 
value of 35.1 instead of 29.6. Some of the discrepancy between curves A 
and C around t = can be ascribed to the assumption of exponential 
service. Better agreement in this region could be obtained by taking the 
model to have the (almost) true service distribution (20). The proce- 
dure for doing this is indicated by (19) but the task of carrying through the 
work seems to be difficult. Again, it should be possible to use the "Erlang 
service" approximation 4b 2 t exp (-2bt) for B'(t) in (19) and also to 
solve the corresponding first passage problem [i.e., solve the equations 
corresponding to (8) and (12)]. However, this was not attempted. 

Some idea of the change produced in G(t) when exponential service is 
replaced by other kinds of service may be obtained from Fig. 2. The 
curves of Fig. 2 show G(t) for the simplified model based upon the re- 
sults of Appendix B. Exponential service (curve D) and constant service 
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time (curve E) are assumed, and Poisson arrivals are taken for both 
cases. The average arrival rate is a = k/T = 18/6.35 = 2.84 per hour 
and the average service time is /;"' = 1.46 hours. Substitution of these 
values in the expressions for G(t) given in examples (a) and (b) of 
Appendix B give curves D and E, respectively. 

The expanded scale at the top of Fig. 2 shows the behavior of G(t) 
around / = 0. Both distributions predict the same average innage 
length, namely 



h = [e ab - \}a" 1 = 22.0 hours 



(23) 



which is (43) in Appendix B. The discrepancy between 22.0 and the 
value ti = 29.6 given by the model of Section III shows the shortcomings 
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Fig. 1 — Innage length distribution for 18-satellite system with random polar 
orbits (Maine-Western Europe). 
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Fig- 2 — Innage length distributions. 

of the simplified model. Nevertheless, it appears that the difference in 
shape between curves D and E illustrates the change in G(t) produced 
by the different kinds of service. Support for this belief comes from the 
fact that curves A and D, both of which correspond to exponential 
service, have the same shape. 

In view of the inaccuracies of the models and of the relatively good 
agreement shown by the curves B and C of Fig. 1, it seems that the 
simple exponential approximation for G(t) is quite good. 
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APPENDIX A 

Probabilities Associated with Alternating Sequences 

Expressions are recorded here for the Laplace transforms of two eon- 
ditional probabilities. These transforms are of use in establishing Takacs' 
result (19). 

Consider a sequence comprised of two kinds of intervals which alter- 
nate with each other (for example, innages and outages). Let p a (t), 
Pb(t) be the probability densities for the lengths of the two types of 
intervals, and let their respective Laplace transforms be a(s), 0(s). 
Suppose that the interval lengths are independent, let an interval of 
type (a) start at time 0, and let P aa (t) be the probability that an interval 
of type (a) is in progress at time t. Then the Laplace transform of 
Pa„{t) is 



/; 



l Pa*(t)dt= r /~;^> (24) 

s[l - a(s)p(s)\ 

Similarly, if an interval of type (b) starts at time 0, the Laplace 
transform of the probability Pi,„(t) that an interval of type (a) is in 
progress at time / is 

[1 - tt(«)lfl(*) , . 

s[\ -a(s)(3(s)}' K J 

These results are reminiscent of the relations between generating func- 
tions in the theory of recurrent events (Ref. 4 ; ch. 12). 

The expression (24) may be obtained by noting that 

P„„(0 = l»r(/, (1 > t) + PrC-i + hi + t a2 > t 

and /.„, + hi < t) + • • • 

= PrU,, > + [l»r(/,„ + /,„ + /„ 2 > 

- Pr(/,„ + ki >/)]+••• 

where l ni , / h i , /,,« , • • • are the lengths of the successive (a) and (b) 
type intervals. The Laplace transforms of Pr(/„i > t), Pr($„i + hi > 
/), • • • are [1 — a(s)]/s, [1 — a(s) /3(s)]/s, • • • , and summation gives 
(24) when | a(s) /3(s) | < 1 (as it certainly is when Rc(s) > 0). Expres- 
sion (2o) may be obtained in a similar fashion or by convoluting Pb(t) 
and P aa (t). 

When p„(t) = a exp (— at), the condition that a type (a) interval 
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start at time may be replaced by the condition that a type (a) interval 
be in progress at time 0. In this case P aa (t) is also the probability that 
a type (a) interval is in progress at time t, given that one was in progress 
at time 0. From (24) and a(s) = a/ (a -f- s), the Laplace transform of 
this probability is 

1 (26) 



s + a — a(j(s) ' 



APPENDIX B 



Sijstems with an Infinite Number of Servers and Recurrent Inputs 

This appendix will be concerned with systems containing an infinite 
number of servers and a "recurrent" input, i.e., an input in which the 
lengths of the intervals between successive arrivals are independent of 
each other and of the state of the servers. In many respects these sys- 
tems are simpler than the limited source input introduced in Section 
III. Although they do not represent the satellite system as well, their 
greater simplicity enables us to estimate the shape of the innage dis- 
tribution for cases which are difficult to handle by Takacs' result (19). 

In the following list of results, A(t) is the distribution function for 
the distances between the arrivals of the recurrent input and B(t) is 
the service time distribution for each one of the infinite number of 
servers. The expected distance between arrivals is a -1 and the expected 
service time is b~ . 

b.i The Conditional Probability Pio (t) 

Let all of the servers be idle at time -0 and let the first customer 
arrive at time making one server busy at +0. Denote by Pi„(t) the 
conditional probability that n servers are busy at time t. Consideration 
of the first arrival following time +0 leads to an integral equation 
which (in theory) may, be solved for Pio(0> namely 

Pl0 (t) = [1 - A(l)]B(t) + B(t) I A'(t - v)PMdv, (27) 

where A'(u) = dA(u)/du. A corresponding equation for the generating 
function for P in (t) is given by (44). 

Example (a). For Poisson input 1 - A(t) is exp (- at). This cor- 
responds to an unlimited source input. Substituting in (27), multiplying 
through by [exp (a/)] B(t), and differentiating with respect to t gives a 
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differential equation for Pu(t). Using the fact that P 10 (0 — » B(0 as 
t — > to fix the constant of integration leads to 



P,o(0 = B(0 exp 



-a/ fl [1 -B(r)]drl, 



(28) 



a result given by Refs. 5 and 6. The work may be simplified by starting 
with the assumption that Pio(0 is of the form B(t)P(t). 

Example (6). For regularly spaced arrivals, A' (a) = 8(u — a~ l ) where 
8(t) is the unit impulse. Equation (27) then gives 



p 10 (0 = B(t), o < t < <r' 

P 10 (/) = B(t)P 10 (t - a' 1 ) t > a" 1 

P,o(0 = B(t)B(t - a~ x ) a~ l < t < 2a" 1 

P l0 (t) = B(t)B(t - a'^Bit - 2a _1 ) 2a" 1 < t < 3a" 1 



(29) 



and so on. 

Example (c). The case when A(t) is arbitrary and 1 — B(t) = 
exp (— t) has been studied by Takacs 7 (see also Ref. 3, p. 33 et seq.). 
Multiplying (27) by exp (— st), integrating t from to °o, and intro- 
ducing the Laplace transforms a(s), 6(s) of A'(t), Pi a (t) leads to a re- 
currence relation between 0(s) and 0(s + 1) which in turn gives 

0(8) = ~ ~ , , *, r + 



8 (s+l)(l-o.) (s + 2)(l- a .)(l- fti+1 ) 

(30) 

<x*+ias+2 , 

(S + 3)(1 - «.)(1 - «H-l)(l - «.+2) ' 

where a B+ „ is written for a(s + n). When 0(s) is known Pu(0 may be 
obtained by inversion. 

b.2 77ic Pwsj/ Period Distribution G(t) for Poisson Arrivals 

Let —G'(t) be the probability density for the lengths of the busy 
periods (corresponding to innages) and consider the case of Poisson 
input and arbitrary service. The Laplace transform y(s) of —G'(t) is 
given by 

_ (s + a)d(s) 

yW 1+00(8) {M) 

where 8(s) is the Laplace transform of Pio(t). 
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My original derivation of (31) was based on taking the Laplace trans- 
form of 



_(?'(/) -/(*) - a [ dsP w (t- s)f(i 



+ a 2 f dsP l0 (t - s) f drl\(s - r)f(r) - 



(32) 



'o *'o 

/(*) = Pio'(0 + flPia(0, (83) 

where P w (t) is given by (28) and /(0 dt is the probability that the 
system will jump from state 1 to state in (*, t + dt), given an arrival 
at time which ends an idle period. The series (32) was obtained by 
an application of the method of inclusion and exclusion. Subsequently, 
Takacs obtained a formula equivalent to (31) by an elegant method based 
upon results of the type stated in Appendix A. When this method is 
applied to obtain (31) the type (a) intervals are taken to be idle periods 
(outages) with p a (t) = a exp (- at), a(s) = a/ (a + s). The type (b) 
intervals become the busy periods so that pb(t), 0(s) become —G'(t), 
y(s), respectively; and Pba(t) goes into P i0 (t). Expression (25) of 
Appendix A then says that the Laplace transform of P w (0 is 

,(.) = T (s) (34) 

s -f a — ay{s) 

from which (31) follows. 

Example (a). For Poisson arrivals and exponential service with 
B = 1 - c~ b ', (28) becomes (Ref. 3, p. 26) 

P 10 (i) = (1 - e- bl ) exp [- p + p c~'"\, P = atr\ 

The change of variable y = «•"*' carries the integral for the Laplace 
transform of P w (t) into 

bd(s) = F(z) - F(z+ 1), z = 8b' 
where 



i-i 



F{z) =e- p fif-Vdy 



1 P _ I P _ . . . (35) 

z z(z + 1) "•* z(z+ l)(z + 2) 



= 2 - [1 _ pF(2 + i)] = E-^)- 
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The Laplace transform of G(t) may be shown to be 

1 - 7(g) = F(z + 1) = 1 __ 1 
s sF(z) azF(z) a 

and inversion gives 



(36) 



G(t) = E *'"" , p = a6- J (37) 

where F'(2) = dF{z)/dz. The zeros of F(z) are real and (a) occur at 
Zq , z\ , 2 2 , • • • , (b) are such that, — 1 < z a < 0, — 2 < z x < — 1 , • • • , 
etc. and hence lie between the poles of F(z) at 0, —1, — 2, • • • , (c) may 
be computed by successive approximations with the help of the last 
series for F(z) in (35). 

A power scries for G(t) may be obtained by expanding [1 — y(s)]/s 
in powers of 1/s and then replacing s~" _1 by t n /n\. The same series may 
be obtained from the corresponding series (5.3) for the limited source 
case by letting k — > °o , a — > in such a way as to keep ka equal to a. 
Replacing a by pb then gives 

ril , , bt ( P + \)(bt)* (p+±p+ \){btf 



2! 3! 

(p' + y + ii P + i)(6/ 4 ) 

4! 

Example (b). For Poisson arrivals and constant service time, B(t) 
jumps from to 1 at time / = & _ . Equation (28) shows that Pio(0 
jumps from to exp ( — ab~ ) at t = b~ . The Laplace transforms are 
readily computed and the one for (l(t) gives 

(!U)= ^i - ■ -^ d *> C> °- (39) 

Taking c large enough, say c = b, to make \x\ < 1 where x is as* 1 
exp [— (a • + 8)6" ] and expanding the denominator in powers of x leads 
to 



G(«) =! + £(- 



pc 



(6/ - ;0" _1 , (W - 



_ (n - l)!p 



-n)" ~\ 
nl J' 



+ VW _.," / . (40) 



Here, as in example (a), p = ab ' and the upper limit of summation is 
determined from N < bt < N + 1. In particular, 



(41) 



^ e 
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0(0) - 1, Gib' 1 - 0) = 1, Gib' 1 + 0) = 1 - e'" 
G{2b~ l ) = 1 - (1 + p) e~ p 

G(3b~ l ) = 1 - (1 + 2p)e~ p + pU + |) 

For large values of t 

o(0 ~ ( ;; ± p \ e - w (42) 

p(l + (To) 

where <r is the rightmost root of 

a + p e - """ = 0. 

When p is large a is approximately — p e _p . 

The innages corresponding to examples (a) and (b) have the same 
average length, namely 

i- = [ e ab ~ l - l]/o. (43) 

To see this, note that from traffic theory, or by letting t — > oo in (28), 
the fraction of idle time is p = 1 - q = exp (— ab' 1 ). The average 
outage time is t Q = I /a, and (43) follows upon using the relation (1) 
between U ,t , and q. 

b.3 Miscellaneous Results for Recurrent Input 

Except for the case in which arrivals occur at multiples of some fixed 
spacing, P w (t) approaches the steady state probability p as t — > », 
and the Laplace transform 6(s) of P w (t) has a simple pole at s = 
with residue p • 

If Pio(t) tends to a periodic function the residue gives its average 
value. Application of this result to the case of regularly spaced arrivals 
and exponential service with B(t) = 1 — e~' leads to the rather curious 
expansion 



L 



i °° 
E 

N=0 



II (1 -.r- T -)dr 

n=0 

= 1 "hTT-L 1 ~2(.r- 1) + 3(;r- l)(x« - 1) J' 

Both sides represent the average value of the periodic function to which 
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Pio(t) tends as t — > oc . The integral on the left (with x = exp a~\ x > 1 ) 
is the average value as computed from (29), while the series on the right 
is the residue of 0(s) at s = obtained by setting a(s) = exp (- saT 1 ) 
in (30) and letting s — > 0. 
The generating function 

;i=0 

for the conditional probability P in (t) that state n exists at time £, given 
that an arrival at time ends an idle period, [see (27)] satisfies the 
integral equation 

Pi(*,0 = [•<• + (! -x)B(t)] 

(44) 



1 - A(t) + / o .4'(/ - w)P,(a:,») < fol. 



A formal step-by-step solution may be obtained by introducing the 
binomial moments M n (t) defined by 

Pi(*,0 - E(a- l)"3/.(0, 



sO 



M n «) = E " )Pu(0. 

The value of il/ () (<) is one and the higher-order moments satisfy integral 
equations obtainable from (44). When the Laplace transforms of 
A'(t), B'(t), M n (t) are denoted by a(s), (3(s), Mn (s) it is found that 
the integral equations lead to 

s[l - a(s)] 

Hn{s) = — . / — _ dz, n > 1. (45) 

2m Je-ico s — z 1 — a(s) 

The singularities of a(z)n n _ l (z) are supposed to lie to the left and those 
of [1 - (3(s - z)]/(s — z) to the right of the path of integration, s 
being chosen so as to make this possible. In theory, the successive values 
of y. n {s) may be obtained step by step and thus ultimately lead to an 
expression for P\(x, t). For exponential service the integrals in (45) 
may be evaluated and lead to results given by Takacs 7 (see also Ref. 3, 
p. 33). 
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